THE CONTACT GEOMETRY OF THE RESTRICTED 3-BODY PROBLEM 
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Abstract. Wc show that the planar circular restricted three body problem is of restricted contact 
type for all energies below the first critical value (action of the first Lagrange point) and for energies 
slightly above it. This opens up the possibility of using the technology of Contact Topology to 
understand this particular dynamical system. 



1. Introduction 

Starting with the fundamental contributions of Gromov and Floer, holomorphic curve techniques 
had a great impact on our understanding of the dynamics of Hamiltonian systems. The three body 
problem is a Hamiltonian system with a particular intriguing dynamics. During the visit of the 
first two authors to the IAS in Princeton, Helmut Hofer raised the question of whether holomorphic 
curve techniques can be applied to the three body problem. 

The power of holomorphic curve techniques relies on a particular compactness result, called 
Gromov compactness. In order to prove Gromov compactness the Hamiltonian system has to meet an 
additional convexity assumption. Namely a Liouville vector field transverse to energy hypersurfaces 
has to exist. We prove in this paper that this assumption holds true for the planar circular restricted 
three body problem for energy values below and slightly above the first critical value. 

The restricted three body problem assumes that two massive primaries, which we refer to as 
the earth and the moon, rotate around each other according to Kepler's law on ellipses. We are 
interested in the dynamics of a third massless object called the satellite. Since the satellite is massless 
it does not influence the propagation of the primaries, however the primaries attract the satellite 
according to Newton's law of gravitation. If the ellipses are circles then the restricted three body 
problem in rotating coordinates admits an integral, which we abbreviate by H. This special case of 
the restricted three body problem is called the circular restricted three body problem. If in addition 
the satellite stays on the plane spanned by the earth and the moon the problem is referred to as the 
planar circular restricted three body problem. 

Energy hypersurfaces for the planar circular restricted three body problems are three dimensional. 
However, they arc not compact. This is due to collisions of the satellite with one of the primaries. 
However, it is well known in celestial mechanics that two body collisions can be regularized. It was 
observed by Moser that the regularized Kepler problem coincides with the geodesic flow of the two 
sphere. In particular, energy hypersurfaces are diffcomorphic to MP 3 . In Moser's regularization the 
position and momenta variables are interchanged and the point at infinity on the sphere corresponds 
to collision orbits whose velocity explodes. 

For low energy levels c C M the energy hypersurface S c = _ff _1 (c) for the planar circular restricted 
three body problem has three connected components. Indeed, the satellite either has to stay close 
to the earth, close to the moon, or be far away. We denote by Sf and the connected component 
close to the earth respectively close to the moon. Via Moser's regularization these components can 
be compactified again to compact manifolds E c and E c which are both diffcomorphic to RP 3 . 

On the axis between earth and moon there is a critical point of the energy called the first La- 
grange point, which we denote by L\. If the energy level c becomes higher than the action of the 
first Lagrange point H{L{), then the satellite is in principle able to cross from the region around 
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the earth to the region around the moon. The energy hypersurface now consists of only two con- 
nected components, one bounded which contains the regions around the earth and the moon and 
an unbounded one. We abbreviate by E^' M the bounded component. Topologically this bounded 
component is just the connected sum of the two bounded components we had before. In particular, 
its regularization E c ' is diffeomorphic to RP 3 f/=RP 3 . 

To state our main theorem we need the following convention. 

Definition 1.1. Let E be a 2n+ 1-dimensional manifold. A Hamiltonian structure on E is a closed 
two-form whose kernel defines a one-dimensional distribution. A contact form on E is a 1-form A 
such that A A (d\) n is a volume form. If u is a Hamiltonian structure on E, then a contact form A 
is called compatible with u if dX = uj. 

Remark 1.2. Abraham and Marsden use in [5] different conventions. A Hamiltonian structure is 
in their language a contact form, where a contact form in our sense corresponds to an exact contact 
form. 

Moser's regularization is actually symplectic and therefore the regularized hypersurfaces are nat- 
urally endowed with a Hamiltonian structure. The main result of this paper is: 

Theorem A. For c < H(L\) both connected components E^ and E^ admit a compatible con- 
tact form A. Moreover, there exists e > such that if c € (if(Li), H(L\) + e) the same assertion 
holds true for E,f ' M . 

Remark 1.3. The drawback of Theorem A is that we only can prove the contact condition until 
slightly above the first critical value. There is no particular reason that the contact condition should 
fail for higher energy values and we actually expect that for every energy level the hypersurfaces are 
contact. However, this is an issue for further research. 

Remark 1.4. Moser regularization actually endows the compact components of the energy hyper- 
surfaces with a strong symplectic filling. In particular, the contact structures obtained in Theorem A 
are of restricted contact type and tight. 

Since there is only one tight contact structure on RP 3 as well as on RP 3 #RP 3 we obtain the 
following corollary of Theorem A. 

Corollary 1.5. For c < H(L\) the contact structures (E c ,kerA) and (E c , kerA) coincide with 

the tight RP 3 and for c G (if (Zi), H{L\) + e) the contact structure (E c ' , ker A) coincides with the 
tight RP 3 #RP 3 . 

We note that in general it is a very difficult problem to decide if a given Hamiltonian structure 
admits a compatible contact form. A notorious concrete case is that of Hamiltonian structures arising 
from Hamiltonians which contain "magnetic terms", i.e., terms which are linear in the momentum 
variables (cf. [HIE])- This is precisely the case of the planar circular restricted three body problem 
which picks up a term of the form p\q2 — Vi<l\ due to the rotating frame. 

Before we mention two applications of Theorem A we point out that the strong symplectic filling 
of the regularized energy hypersurface E c is T^ :i S 2 , the unit co-disk bundle of S 2 , and of E c ' is 
the boundary connected sum T^ 1 S 2 \\T^ 1 S' 2 , see Remark 7.4 

The first application concerns existence of periodic orbits. It follows from the fact that the 
contact structures are tight, and Rabinowitz's result on existence of closed orbits on starshaped 
hypersurfaces [17]. The details of this derivation are explained by Hofer in [13] together with the 
fact that the lift of the tight contact structure on RP 3 to S 3 is the unique tight contact structure 
on S 3 . 

Corollary 1.6. For any value c < H{L\), the regularized planar circular restricted three body 
problem has a closed orbit with energy c. 
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Remark 1.7. The assertion of the previous corollary is trivially true for larger values of c since the 
standard Weinstein model for the connected sum contains a periodic orbit. 

The second application concerns leaf- wise intersection points. This notion is due to Moser [TB]. We 
recall that if S C (W, ui) is a regular hypersurface then E is foliated by the leaves of the characteristic 
line bundle kerw|s. We denote by L x the leaf through ieS and furthermore by Ham c (W / ) the set 
of Hamiltonian diffcomorphisms generated by a Hamiltonian function with compact support. Then, 
by definition, x € S is a leaf- wise intersection point with respect to tp € Ham c (W) if ip(x) € L x . 

Corollary 1.8. Let e be the same as in Theorem A. Then for c < H(L\) + e, c =/= H{L{) and any 

tp e Ham c (T*5 2 ) the regularized energy hypersurface E c ' carries a leaf-wise intersection point 
with respect to ip . 

Proof. Since the contact structures are of restricted contact type we can apply the theorem by 
Cieliebak-Frauenfelder-Oancea [B] and obtain that the Rabinowitz Floer homology of the energy 
hypersurface as defined by Cieliebak-Frauenfelder [5] is isomorphic to the direct sum of symplectic 
homology and cohomology (except possibly in finitely many degrees.) By Cieliebak's [4] computation 
of the symplectic homology of a boundary-connected sum, see also McLean [TJJ Theorem 2.20] for 
an improvement concerning the ring structure, we know 

SH^T^S^T^S 2 ) S SH^T^S 2 ) x SH^T^S 2 ) . 

Combining this with the work of Viterbo [19] . Salamon- Weber [18] . Abbondandolo-Schwarz [TJ 

SH^T^S 2 ) S H^AS 2 ) 

where AS 2 denotes the free loop space of S 2 we see that in both cases, c < H(Li) and c < H(Li) + e, 
Rabinowitz Floer homology of the energy hypersurface is infinite dimensional, in particular, non- 
zero. Thus, work by Albers-Frauenfelder j3j Theorem C] implies the Corollary. □ 

Remark 1.9. We remark that Theorem A and all its corollaries are valid for any mass ratio between 
the two primaries. 

2. Elements of 3-dimensional contact topology 

In the proof of Theorem A it is convenient, although not necessary, to use some contact topology. 
In doing so, we can regard the dynamical problem from another perspective. Let us begin by defining 
some key notions. 

Definition 2.1. A contacts-manifold £) is called overtwisted if(M,£) admits an overtwisted 
disk, i.e., an embedded disk A such that £ p = T p A for all p € OA. Contact 3-manifolds that are not 
overtwisted are called tight. 

The contact topological behavior of tight and overtwisted contact manifolds differs very much. 
Roughly speaking, one could say that tight 3-manifolds are much more rigid than overtwisted ones. 
Finding tight contact structures is a hard problem in general, but there is a simple criterion due to 
Eliashberg and Gromov, see [TT] and |12) . 

Theorem 2.2 (Eliashberg, Gromov). Suppose (M, is a symplectically fillable 3-manifold, then 
(M, £) is tight. 

By symplectic fillability we mean the following. 

Definition 2.3. A symplectic manifold (W,u>) is called a strong symplectic filling for (M,£) if 
dW = M and there is a Liouville vector field X pointing outward along M such that £ = kei(ixuj) \m- 

Some uniqueness results are known about tight 3-manifolds. Note that the general situation can 
be much more complicated, but the following suffices for our needs. 
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Theorem 2.4 (Eliashberg) . The closed 3-manifold M.P 3 admits a unique tight contact structure up 
to isotopy. 

This unique tight contact structure is, in fact, symplectically tillable; its filling is T*S 2 . Alterna- 
tively, we can see (RP 3 ,£o) as the contact manifold obtained from (S 3 ,^o) by dividing out by the 
antipodal map. 

Finally, let us point out that tight manifolds admit a unique prime decomposition, see [10] . 

Theorem 2.5 (Ding-Geiges) . Every non-trivial tight contact 3-manifold (M, £) is contactomorphic 
to a connected sum 

(M,W#-#(M*,f fc ) 

of finitely many prime tight contact 3-manifolds. The summands are unique up to order and con- 
tactomorphism. 

By prime we mean that a manifold is not decomposable in a non-trivial connected sum; non- 
triviality is needed to exclude connected sums with (S 3 ,£,o)- 

2.1. Weinstein model for connected sums and symplectic handles. Let us briefly recall the 
Weinstein model for connected sums. Consider (R 2n+2 ,u;o = dx A dy + dz A dw), where x,y G K™ 
and z, w € K.. Observe that the Liouville vector field 

X = - (xd s + ydg) + 2zd z - wd w 

is transverse to level sets of / = x 2 + y 2 + z 2 — w 2 provided / ^ 0. This observation is the starting 
point to show that connected sums of contact manifolds carry a contact structure. 

Let (Mi,£i) and (Ma, £2) be contact manifolds, and choose Darboux balls Di C Mj. To perform 
the contact connected sum, observe that we can also embed D1UD2 in the two connected components 
°f !)• hi f ac t, ° nc could regard / _1 ( — 1) as two disjoint Darboux balls embedded as contact 

hypersurfaces in M. 2n+2 . Now connect these two Darboux balls by "interpolating" / _1 (— 1) to 
This is the neck of the connected sum. Note that this interpolation process is possible, since 1) 
and / -1 (1) are close to each other for large w. The Liouville vector field X is transverse to this 
neck, and hence the connected sum is contact. 

Remark 2.6. Alternatively, one could use the Liouville flow to "map" / _1 (— 1) to In this 

process, some isotopy is still necessary, though. 

2.2. Symplectic handles. Suppose we want to perform symplectic handle attachment. To do so, 
consider the symplectic 1-handle, 

SH := {-2 < / < 2}, 

and let W be a (not necessarily connected) symplectic cobordism W with boundary dW — M . Take 
Pi and P2 points in M along which we shall attach the handle. 

By Darboux's theorem we find neighborhoods Di of pi that are (strictly) contactomorphic to 
(R 2n + , ao). That means that we can also find a strict contactomorphism ipi U if 2 from D\ U D2 to 
a compact subset of the two components of the level set / _1 (— 1). 

Furthermore, W is a symplectic cobordism, so a collar neighborhood of the boundary is symplecto- 
morphic to ([— 1, 0] x M, d(e t a)). We can construct a symplectomorphism from a collar neighborhood 
of Di to a compact subset of SH by gluing flow lines of the Liouville vector field on W (dt in a 
collar neighborhood of the boundary) to flow lines of the Liouville vector field on SH. 

In other words, the symplectomorphism is of the form 

fa : [-1, 0] x Di C [-1, 0] x M — > SH 

(t,x) .— > Flf(ipi(x)). 

Remark 2.7. The above formula is not quite accurate, since the amount a flow lines changes the 
value of / in SH depends on the the precise starting point. In other words, ipi maps the collar 
neighborhood [—1,0] x Di into a compact subset of SH, see Figure [l] for a sketch. 
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Define SH C to be the compact subset of SH consisting of point x that lie on a flow line of X 
through ip(0,Di). Define the symplectic manifold obtained by handle attachment by 

W := W\JSH C / ~, 

where we identify (t,x) G [—1,0] x Di c W with p in SH C if and only if p = ip(t,x). See Figure [l] 
for an illustration of the symplectic handle attachment. 




Figure 1. A symplectic handle 



3. MOSER REGULARIZATION 

In this section we recall the regularization introduced by J. Moser in [15] . This regularization 
will be used in subsequent sections, but here as a warm-up problem we apply it to the following 
integrable Hamiltonian on T*(R 2 \ {(0, 0)}): 

(1) H(p,q) = ^\p\ 2 - +piq 2 -p 2 qi- 



This Hamiltonian corresponds precisely to the case when we set fj, = in ( 10 1 and it describes the 
Kepler problem in a rotating frame. 

The first step is to change the independent variable t to s by setting 

dt 



and consider a new Hamiltonian K defined by 

H = j— r + k. 

m 

Here k is the energy level of H that we are interested in regularizing. The observation now is that 
for points at K = (which corresponds to H = k) one has 

H _K p _ Kq 

P - W\ ' 9 ~ k\ ■ 

Then the orbits of H with energy k and time parameter t correspond to orbits in K — with time 
parameter s. We now compute K and we obtain: 

K = \{\p\ 2 + l)\q\ + \q\(piq 2 -pa«i) + (-* - 1/2)| 9 | - f. 
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The second step consists in performing the canonical transformation p = —x and q = y, so that now 
we would like to understand the zero energy level of: 

(2) K = X -{\x\ 2 + l)\y\ + \y\(-x m + x 2Vl ) + (-k - l/2)\y\ - 1. 

The last step is to introduce stereographic projection. Using the same notation as in [T5J, we let 
£ = (£o)£i>£2) be a point in R 3 with norm one, so it represents a point in S 2 . A tangent vector 
r) € T^S 2 is written as r\ = (?7o> ^li %)> with inner product (£, 77) = 0. We shall identify TS 2 with 
T* S 2 using the standard metric on S 2 . The transformation induced by stereographic projection is 
described by the equations: 

(3) Xk = -^-, k = 1,2, 

1 - so 

(4) lfc=»fr(l k = l,2. 

This is a symplectic transformation between T*M. 2 and T* S 2 and we also have the relations: 

(5) (N 2 + i)M _ |y| 

We now apply these formulas to the Hamiltonian K in (J2| and we find after a calculation that our 
system goes over the hypersurface in T* S 2 given by: 

M + (i - £o)M(r7i6 - mb) + (1 - £o)M(-fc - 1/2) = 1. 

To make this smooth at the zero section we consider the Hamiltonian on T*S 2 

(6) Q(f, r?) = \\v\ 2 [l + (1 - &)(-* - 1/2 + - r/ 2 Ci)] 2 

and we are interested in the hypersurface Sfe := Q _1 (l/2). In a moment we will show that the 
hypersurfaces are all of contact type for the relevant range of energies (k < —3/2), but first a 
remark. 

Remark 3.1. If we had started without the term p\q 2 — P2I1, which corresponds to the rotating 
frame, the Hamiltonian Q would be 

Q(e^) = ^M 2 [i + (i-e )(-fc-i/2)] 2 . 

For k = —1/2, which is the energy considered by Moser in [15] . the system is precisely the geodesic 
flow of the standard metric on S* 2 . For k < (the relevant range of energies for the Kepler problem) 
we also obtain a Riemannian metric on S 2 given by 

1 

9k ~ [l + U-£o)(-fc-l/2)] 2 9st 

where g s t is the standard metric on S 2 . An elementary calculation (that we omit) shows that the 
curvature of gu is constant an equal to —2k, so up to scaling, g k is isometric to g s t- 

Returning to the Hamiltonian ([I]), we rewrite it as: 

H(q,p) = 1 ((p! + q 2 f + (ft, - qi f) + U(q), 

where U is the effective potential. If we let r = |g| then 

U(r) = -l/r-r 2 /2. 



The boundary of the Hill's region (see ( 12 ) for the definition) with energy k is given by 

fc + l/r + r 2 /2 = 

so we look at the cubic 



r 



3 - L 2kr + 2 = 0. 
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It follows that there is a compact Hill's region iff k < —3/2. So for k < —3/2 the compact Hill's 
region is given by r < tq where ro(fc) is the smallest positive root of the cubic (ro(— 3/2) = 1). Thus 
k < —3/2 is the range of energies for which has a compact connected component containing the 
zero section. We now show: 

Proposition 3.2. is starshaped for k < —3/2. 

Proof. We shall check that the Liouville vector field X = rjd v is transverse to . Since everything is 
rotationally symmetric about the 6~ ax i s j it is enough to consider the case of 6 = and £q +£§ = 1- 
So we look at Efc over these points only: 

(7) |»?|[l + (l-£o)(a + »?i6)]=l 

where a := —1/2 — k. Now for A near 1 we let 

/(A) :=A|r ? |[l + (l-£o)(a + A7 7 i6). 

We wish to check that /'(l) ^ for points satisfying Q. After a calculation we find 

l + (l-eo)(a + 27 ?1 6) 
J[ > l + (l-€o)(a + »?i&) ' 

Thus /'(l) =0 iff 

,_ a _ l + (l-g )a 

77142 2 2(1 -6) 2(1-6) ' 

Note that at the point with 6 = 1 there is nothing to check and without loss of generality we may 
assume 6 > 0. Using ^ we see that still assuming /'(l) = 



Clearly \r]\ > \r]i\, thus 

(i + (i-6)q) 2 
46(i-6) 

Now, if we let 6 = cos $ an d 6 = sm $ we can rewrite the above as 



< 1. 



(l + (l-cosg) a ) 2 
5a(0) - 4sin#(l-cos0) " L 



Claim. If a > 1, then g a > 1. 
We introduce the function 



, , (2-0 4 (i-2)" 1 

ft: (-1,1) — > R, 1 1— > 



16(l-t 2 )(l-t) 2 16(t + l)(t-l) 3 ' 



Note that for 8 E (0, ir) we have 

(8) ft(cosfl) =gi(0) 2 . 

The derivative of ft is given at t € (—1,1) by the formula 

3i(i-2) 3 



ft'(i) 



8(i + l) 2 (t-l) 4 - 

In particular, the only zero of ft' in (—1,1) is located at t = 0. Because liiiit_»±i ft(£) = oo we 
conclude that ft attains its global minimum at t = 0. Since ft(0) = 1 we deduce 

(9) ft > 1. 

Combining ([8| and ([9]) we conclude that 

9i > 1. 

Since for a > 1 it holds that g a > <?i the assertion of the Claim follows. 
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It follows that for k < —3/2, /'(l) ^ and we have transversality of the Liouville vector field 
with Efc. 

□ 

4. The circular planar restricted 3-body problem 

Assume that we have two massive bodies called the primaries, which we will denote as Earth and 
Moon, which propagate according to Newton's law of gravitation on circles around their common 
center of mass. If wie is the mass of the earth and tum is the mass of the moon we denote by 

mM c rnii 
/! = ■ G 0, 1 

m E + m M 

the proportion of the mass of the moon on the total mass. In the following we normalize our problem 
such that the total mass satisfies 

m E + m M = 1- 

In an inertial plane spanned by the earth and the moon the position of the earth becomes 

E(t) — (^cost, ^sini) 

and the position of the moon 

M(t) = ( - (1 -M)cosf,-(l -/z)sint). 

Suppose further that S is a third massless body moving in the plane spanned by the earth and moon 
on which the gravitational force of the two primaries act. The body S is referred to as the satellite. 
Since the satellite is massless it does not influence the two primaries. For t g R the time-dependent 
Hamiltonian 

E\ : (R 2 \ {E(t),M(t)}) xM 2 4l 
for the satellite is given at (q,p) e (R 2 \ {E(t),M(t)}) x R 2 by 

rrif s 1, ,2 1-M i 1 

= -h=m\-\ q - M (t)\' 

i.e. the sum of the kinetic and potential energy of the satellite. Here the superscript i indicates 
that we consider an inertial coordinate system. Instead of considering an inertial coordinate system 
we can alternatively also study the propagation of the satellite in moving or synodical coordinates 
where the earth and moon are at rest. We choose the coordinates in such way that the position of 
the earth and moon satisfy 

£=( M) 0), M= (-(1- M ),0). 
The Hamiltonian in the synodical coordinate system for the satellite is given by 

H : (R 2 \ {E, M}) xl 2 4l 

satisfying 

(10) H(q,p) = hp\ 2 - ^ — ^— +piq 2 -Pzqx- 

2 \q-E\ \q-M\ 

The transition from the inertial coordinate system to the moving coordinate system requires a time- 
dependent transformation and we refer to [2] for a derivation of it. It is an amazing fact that in 
the rotating coordinate system the Hamiltonian becomes autonomous, i.e. independent of time. In 
particular, the Hamiltonian H is preserved under the flow. This observation goes back to Jacobi 
and the integral — 2H is traditionally called the Jacobi integral. 

We further note that the Hamiltonian in the synodical coordinate system is not anymore of the 
simple form kinetic plus potential energy but contains the third term piq 2 — Piq\ which is the 
moment map of the 5 1 -action on the cotangent bundle of R 2 induced from rotations of the base. 
One might think of this term as "rotational energy" governing the centrifugal and Coriolis forces 
of the rotating system. In the case where the Hamiltonian is just the sum of kinetic and potential 
energy it is well-known that energy hypersurfaces are of contact type. However, because of the 



CONTACTING THE MOON 



<> 



rotational term the question of whether energy hypersufaces of the Hamiltonian H are contact is far 
from obvious. Indeed, formally the rotational forces are nothing but the Lorentz force of a magnetic 
field (after switching to the effective potential, see below) and we refer to [7J H] for examples of 
energy hypersurfaces which are not contact in the presence of a magnetic field. 

If we complete the squares in the formula for the Hamiltonian H we end up with the expression 



H(q,p) = \([pi+q2? + {p2-qi?) - y 



M I* l |g|2 



The function 
given by 



q-E\ \q-M\ 2 
U: R 2 \ {E, M } -> K 

U(q) = --, — { - | ryy - -\q\ 

\q-E\ \q-M\ 2 



is called the effective potential. So we obtain the following condensed formula for the Hamiltonian 

(11) H(q,p) = l - (V + q 2 ) 2 + (pa - qx) 2 ) + U(q). 

This expression of the Hamiltonian shows clearly the presence of an exact magnetic field with 
primitive given by the 1-form q 2 dqi — q\dq 2 . 

The effective potential U has five critical points (I 1 , . . . ,£ 5 ) called Lagrangian points. The first 
three critical points, which were already discovered by Euler, are called collinear Lagrangian points 
since they lie on the axis spanned by earth and moon. The first critical point I lies in between 
earth and moon, £ 2 on the opposite side of the moon and £ 3 on the opposite side of the earth. The 
three collinear Lagrange points are saddle points of the effective potential. The two maxima of U lie 
at £ 4 and t which are called the triangular Lagrange points since both of them span together with 
the two primaries an equilateral triangle. Note that if R: R 2 — > M 2 is the reflection along the first 
axis then R(£ 4 ) — £ 5 which is obvious from the fact that the effective potential is invariant under R. 
We remark that the triangular Lagrangian points are also sometimes referred to as Trojan points 
since in the Sun- Jupiter system they correspond to the locations of the Trojan asteroids. 

There is a one-to-one correspondence between critical points of the effective potential and critical 
points of the Hamiltonian. Namely for j S {1, . . . , 5} abbreviate 

L, = {£[,£{,-£{,£{) eR\ 

Then L\, . . . , L5 are the five critical points of H. In case where the mass of the moon is less than 
the mass of the earth, i.e. fi < | the critical points arc ordered by action in the following way 

H(L X ) < H(L 2 ) < H{L 3 ) < H{U) = H(L 5 ). 

If c € K we abbreviate the energy hypersurface for the Hamiltonian H at energy level c by 

s c = /r 1 (c). 

We further abbreviate by 

7T : (M 2 \ {E, M}) xM 2 4l 2 \ {E, M] 
the projection along R 2 . The Hill's region is defined as 

(12) K, c = 7r(S c ) cl 2 \ {E, M }. 

This is the region in position space where the satellite of energy c is allowed to stay. Note that since 



the first two terms in ( 11 ) can attain in each fiber arbitrary nonnegative values, we can alternatively 



define the Hill's region as 

(13) K. c = {qeR 2 \{E,M} :U(q) <c}. 
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If c < H(Li), then JC C contains three connected components, an unbounded one and two bounded 
ones in whose closure lies the earth respectively the moon. We denote the bounded components by 
JCf and Kf such that 

E g cl(/Cf ), M e cl(ACf ). 

We further abbreviate by 

sf = 7 r- 1 (/cf)ns c , sf = 7r - 1 (/cf)nE c 

the corresponding connected components of the energy hypersurface S c . 

5. Below the first critical level 

Assume that c < H{L\) and recall that C i? _1 (c) is the connected component of the energy 
hypersurface around the moon. We introduce the Liouville vector field 

(14) X = {q-M)^. 

The main result of this section is 

Proposition 5.1. X intersects transver sally. 

To prove the proposition we need three Lemmas. The Lemmas need several calculations the details 
of which arc relegated to the Appendix. To formulate the first Lemma we make the convention that 
(p, 9) are lunar polar coordinates, i.e. polar coordinates centered at the moon, such that the effective 
potential becomes 

U ^ = --n~ I 2 a \ V +P C ^ 1 A U 1 ^) 2 - 

P \J p A — 2/jcos + 1 ^ * 

For p e (0, 1) we define 

U p = U(p,-)€C°°(S 1 ,R). 
The first Lemma we need is the following assertion. 

Lemma 5.2. For p e (0, 1) the function U p attains its minimum at 6 = 0. 
Proof. The derivative of U p is given by 

U'{6) = (1 - p)psme(- 1 -3- - 1 



(p 2 - 2pcos8 + l)' 

We see that 9 = and 9 = tt are critical points of U p . The second derivative at these points is given 

by 

^(0) = (1-M)(( T ^)3- 1 )> 



respectively 



W = -(l-M)(^-l)>0. 



We conclude that U p attains a local minimum at and tt. The remaining critical points of U p 
necessarily satisfy the condition cos 9 = % ■ In the two arcs bounded by the minima there is precisely 
one such critical point, hence it has to be a maximum. Therefore the only local minima lie at 9 = 
and 8 = n. We compare the values of U p at them and obtain 

2( 1 - P)P 3 
P)(l + P) 

We conclude that U p attains its global minimum at 9 = 0. □ 



U p (0)-U p (tt) = -—^ — -<(). 
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Abbreviate 

d := \f - M\ 

the distance from the moon to the first Lagrangian point, i.e. the critical point of the effective 
potential which lies on the earth-moon axis between the two primaries. We further introduce 

B = {q G M 2 : \q - M\ < d} 

the ball centered at the moon of radius d. In particular, l x lies on the boundary of B. As a corollary 
of the previous Lemma we obtain the following bound on the component of the Hill's region around 
the moon. We denote by int the interior of a set. 

Corollary 5.3. K,f C int(fl). 

Proof. We first show that 

(15) JC c C\dB = %. 



To see that pick (d,9) G dB. We estimate using Lemma 5.2 



(16) U(d, 9) = U d {9) > U d (0) = U(f) = H (L 1 ) > c. 



By the description of the Hill's region given in ( 13 ) we see that (p, 9) does not lie in JC C and therefore 
( 15 1 follows. We deduce that both sets K.™ n B and n B c are open and closed. Since M lies in 
the closure of /C* 7 , the set IC^f H B is nonempty. Hence, since IC^ 1 is connected, we conclude that 
Kf n B c = 0. This finishes the proof of the Corollary. □ 



The next Lemma we need in order to prove Proposition |5.1| is the following assertion. 
Lemma 5.4. For every q G B - {MJ 1 } it holds that f%) > 0. 
Proof. We prove the Lemma in three steps. 

Step 1: If p G (0,d), then fp(p,0) > 0, i.e. the assertion of the Lemma holds on the earth- 
moon axis between the moon and the first Lagrange point. 

We consider the function u G C°°((0, 1),K.) given by 

u(p) = U(p,0), pG[0,oo). 
Explicitly, the function u is given by 

(17) u(p) = -^- 1 —^- 1 -(p-l + rf. 

p 1 - p 2 

Its second derivative satisfies for p G (0, 1) 

2p 2(1 — p) 

u {p) = j - 7: - K 0. 

r (1 - P) J 

In particular, u is strictly concave. Since the Lagrange point is a critical point of U, we have 
u'{d) = 0, from which the assertion of Step 1 follows. 

Step 2: For p G (0,1) let V p G C^S^M) 6e the function which is given for 9 e S 1 by V p (9) = 
^(p,9). Then V p attains its unique minimum at 9 = 0. 

We first compute the derivative of V p to be 

d ^ V(/O 2 -2pcos0+l) 2 / 
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We first examine the two critical points at 8 — and 8 — n. The second derivatives at these points 
compute to be 

W-a-,)(=^-i)-a-,)(££-i)>o 

and 

w~p-")(=£3?H>° 

so that we conclude that V p has strict local minima at 8 = and 8 = tt. 

We next show that there are no other local minima of V p than 8 = and 8 = tt. For this purpose 
we have a closer look at the function / e C°°([— 1, 1], R) given by 

-2p 2 + + 1 

f(i-) = - tt- 1 ' Te [M- 

(p 2 -2pr + l) 3 

As we just seen above we have 

/(1)>0, /(-1)<0. 

Moreover, the function / is strictly monotone. Hence there exists a unique r € (—1, 1) such that 
f( T o) = 0. Critical points of V p different from 8 = or 8 = tt are precisely points 8 <E S 1 which 
meet the condition cos6> = r . Hence there is precisely one critical point in the interval (0, ir) and 
(tt,2tt). Since V p has local minima at 8 = and 8 — tt wc conclude that these other critical points 
correspond to local maxima. Hence 8 = and 8 = tt are the only local minima. 

To determine the global minimum of V p it suffices now to compare the values at 8 = and 8 = n. 
We obtain 

V p (0) V p (tt) = 2(1 - M) (1 _^~^ p)2 < 0. 
Hence = is the unique global minimum for V p . 



Step 3: Proof of the Lemma. 

By Step 1 it suffices to show the Lemma for 8^0. Hence using Step 1 again together with 
Step 2 we get for p e (0, d] and 8 e S* 1 the estimate 

dU dU 

— (p,8) = V p (9)>V p (0) = —(p,0)>0. 

This finishes the proof of the Lemma. □ 

We continue to denote by B the ball whose radius is the distance from the moon to the first 
Lagrange point and which is centered at the moon. 

Lemma 5.5. For each q e B - {M}, it holds that ^r{q) < — 1. 

Proof. We first compute 

a ^{p,6) = - 2 4 — r(2p 2 - 4pcos0 - 1 + 3cos 2 ^ - 1 

dp 2 ^ ' p 3 (p2_ 2 pcos0 + l)f V ' p ) 

We introduce the function W G C°°{B — {M},R) which is given in lunar polar coordinates by the 
formula 

W(p,8) = ^ f 2p 2 -4pcos(9- 1 + 3cos 2 6»V 

P 3 (p 2 -2pcos0 + l)S V J 

Hence the assertion of the Lemma is equivalent to 

W(q) < 0, V qeB- {M}. 
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We prove this assertion in four Steps. 

Step 1: The distance d from the moon to the first Lagrange point is related to p by the following 
formula 

_ d 5 - 3d 4 + 3d 3 

( j M ~ d i -2d 3 -d 2 + 2d-r 

The distance d is given as the unique critical point in the interval [0, 1] of the function u introduced 
in (17). The derivative of u at p e (0, 1) is 

which gives the relation between p and d as claimed. 

Step 2: For p e (0,1) let W p = W(p,-) E C 00 (S 1 , R) . Let d e S 1 be a point where W p at- 
tains its maximum. Then i? satisfies 

, p 2 - 1 + V-p 4 + p 2 + 1 

cosv = . 

p 

We first compute 



W' p {6) = - 2 3(1 o ^f* ( pcos 2 (9 + 2(1 - p 2 ) cos + p(2p 2 - : 
(p z — 2pcos 8 + 1) 2 \ 

We see that the function W p has critical points at 8 = and 8 = tt. We claim that both of these 
critical points are local minima. We first check this for 8 = 0. The second derivative of W p at 8 = 
is given by 

rtl = 6(l- M )(p 3 -p 2 - P +l) 6(1- M )Q9 + 1) 

_ (1-P) 7 _ (1-P) 5 
This shows that = is a local minimum for W p . Similarly, we get for 8 = tt 

w n, , 6(1- p)(p 3 + p 2 -p- I) 6(1- M )(1 + P) 2 ^ n 

w " {v) = (wF = H-Pf >0 

implying that = 7r is a local minimum as well. Therefore the local maximum has to satisfy 

pcos 2 i? + 2(l- ( o 2 )costf + p(2p 2 -3) =0. 

In particular, 

, p 2 - 1 ± V-p 4 + p 2 + 1 

costf = . 

p 

Since 

P 2 - 1 - V-P 4 + P 2 + 1 < _j 



the assertion of Step 2 follows. 

Step 3: The function W attains its maximum at the boundary of B. 

We argue by contradiction and assume that (po,8q) is a point in the interior of B at which W 
attains its maximum. We first observe that 

(19) = — (po,e ) = -y +3(1- p)(cos8 - po) —„ — -7 ■ 

dp pi (^-2/3ocos0 o + l) 3 

Moreover, at 8q the function W po attains its maximum, so that we obtain from Step 2 the equality 

(20) co S 8 Q = P °- 1 + ^- p ° + p Z + 1 . 

Po 
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In particular, we have 

(21) cos6» -po>0. 



Moreover, taking once more advantage of equality (|20|), we obtain 
3 + 4p cos O - 2/9q - 5 cos 2 O 



2p 4 + 4pg _ io + (10 - 6 P 6V-Po + Po + l 



Pi 

(10 - 6 P 2 ) 2 (- Po + P l + 1) - (2p 4 + Apl - 10) 2 
p 2 (l0 - 2 P 4 _ Apl + (10 - 6p§V-pg + P§ + l) 

20(3 - 8 P 2 + lp\ - 2pg) 
10 - 2p 4 - Apl + (10 - 6p§)V-pg + P§ + l 
40(p o -l) 2 (l-^) 



10 - 2p* - 4p 2 + (10 - 6p§ V-pg + p§ + 1 
> 0. 



Together with (21 1 this contradicts (19). This finishes the proof of Step 3. 
Step 4: Proof that W is nonpositive. 

Since W(q) tends to — oo when q tends to the moon, we conclude that there exists a point (po, 9q) G B 
at which W attains its maximum. By Step 3 we get that 

Pa = d. 

By Step 2 we conclude that 

d 2 - 1 + V-d 4 + <P + l 

COS 0q = . 

a 

Talking also advantage of Step 1 we obtain 



2d 2 (d 2 - 3d + 3) (3 - d 2 - 2 V /= cF + d 2 + l) 

W(Po,0o) = 



(<i 4 - 2d 3 - d 2 + 2d - 1) (3 - d 2 - 2V-d 4 + d 2 + l) 1 d 2 

(d 5 - 2d 4 + d 3 - d 2 + 2d - 1) (6 - 2d 4 - (6 - 2d 2 )V-d 4 + d 2 + l) 

(d 4 - 2d 3 - d 2 + 2d - 1) (3 - d 2 - 2V-d 4 + d 2 + l) i d 2 

The righthand side is now a function which only depends on d E [0, 1]. Plotting this function reveals 
that it is always negative. An analytic argument is included in the appendix. 

□ 

Armed with the previous three Lemmas we are now in position to prove the main result of this 
section. 



Proof of Proposition \5.1\ In lunar polar coordinates the Liouville vector field X is given by 



By ( 11 ) the Hamiltonian in lunar polar coordinates is given by 

H(p, 0) = \((P1+P sin 0) 2 + (pa - p cos + 1 - At ) 2 ) + U{p, 9). 
To prove the proposition we show that 

X{H)\ s m > 0. 
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We first estimate using the Cauchy-Schwarz inequality 



X(H) 



p sin 9(pi + p sin 9) + p cos 9(p 2 — p cos 9 + 1 — p) + p 



dU 
~3p 



> P 



dU 
dp 



\J p 2 cos 2 6 + p 2 sin 2 0\/(pi + p sin 9) 2 + (p 2 — p cos 9 + 1 — p) 2 



Pitt 

= p°~ - p^2{H - U). 
If we restrict this to the energy hypersurface we obtain 

X{H)\ S M>p(^-y/2(c-U) 



The righthand side is independent of the momentum variables so that it is sufficient to check its 
positivity on the Hill's region. Hence we are left with showing 



(22) 



dU 
dp 



y/2(c-U) 



> 0. 



We pick a point {p, 9) € JC^ . By Corollary 5.3 we have p < d. Moreover, by ( 16 ) we have that 



U(d, 9) > c. Since U (p, 9) < c by ( 13 ) we conclude that there exists r £ [0, d — p) such that 

U(p + T,9) = c. 

Hence we estimate 

" T d fdU{p + t,( 



( dU( P ,9) 
\ dp 



> 



f dU(p + T,9) 
\ dp J J dt \ dp 

2 : T du( P + t,9) d 2 u( P + t,9) dt 



dt 



dp 



dp 2 



> 2[ TdU ^ldt 
Ja dp 

= 2(U{ P + T,6)-U(p,9)) 
= 2(c-U(p,9) 



In the second step we have used Lemm a 15. 4| and in the third step we have used again Lemma 5.4 
together with Lemma 5.5 This implies (|22|) and hence the Proposition follows. □ 



Remark 5.6. The above proof actually shows more than the claim of the proposition. For e > 
small enough, the Liouville field is still transverse to level sets with energy E < H(Li) + e away 
from the Lagrange point L\. 

To see this, observe first of all that we can remove a ball around the Lagrange point Bg(Li) to 
divide the level set H = E into a "moon" and an "earth" component. We shall continue to consider 
the "moon" component and argue that the "moon" component is entirely contained in Bd(M). 
Indeed, recall that for fixed p the function U p (9) attains its minimum at U P (Q), so the distance to 
the moon is maximal along the ray 9 = 0. That means that, for e small enough, the maximum 
distance to the moon M occurs for small 9, which we can assume to be less than d provided that e 
and S are small enough. 



With this observation in mind, we can repeat the argument of the proof of Proposition 5.1 to see 
that X is still transverse to the level set H = E for points x that lie in the "moon" component (so 
in particular outside Bs(Li)). 
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6. Regularizing level sets of H and Liouville vector fields 

The goal of this section is to show that we can perform Moser's regularization, and extend the 
Liouville vector field (fl4| to the regularization. We shall use the Hamiltonian 



(23) H(p,q) = l-\p\ 2 +p 1 q 2 -p 2 qi ' 



2"' ~~ \q-q°\ \q-ql[ 

Here q° = (—(1 — fi), 0) is the position of the moon, and q 1 — (fi, 0) the position of the earth in the 
rotating coordinate system. 

Since the only problem in compactness of the components of the level sets corresponding to the 
earth or moon comes from points q with either \q — q \ < e or \q — q x \ < e, we shall only consider 
that situation. Without loss of generality, we shall assume that \q — q°\ < e. 

6.1. Transforming the Hamiltonian. Here we follow the notation introduced in Section [3] In 
order to regularize £, we begin by reparametrizing the flow, i.e. we put 

dt 



\q-q°\ 



and introduce a new Hamiltonian by 



|g-«°l 



q-q° 



q-q 1 



By examining Hamilton's equations we see that the flow of K at energy level corresponds to the 
flow of H at energy level k. The new Hamiltonian is explicitely 

K = \(\p\ 2 + l)\q- <Z°| + \q- q \(pm-p2qi) - (k + \)\q — q°\ — fi — (1 — p) 

Now use the canonical transformation p = —x, q — q° = y. We get 

K = * (M 2 + l)|y| + \y\(x 2yi - x lV2 ) + |y|0r 2<? ° - x iq a 2 ) (k + hy\ (1 - M ) ly J _ - » 

z z \y i q 5i 

Remark 6.1. This transformation reverses the roles of what one would usually expect from the 
position and momentum coordinates; this will become clear in the following. 

Perform the inverse of the stereographic projection given by ([3]), Q and jjjjl. In particular, we 
have 

e > = \y\ = tMi = Jt'U = M^-w- 

We get the new Hamiltonian 

K(Z,v) = M (i + (i - toMm - + (6?? - eq 2 ) ~(k+ i)(l - Co) 

-(1 - M )(l - g ) * j - fi. 

- to) + Zvo + q° - q'\ J 

From this we get the following smooth Hamiltonian on T* S 2 (shifting and squaring the previous 
Hamiltonian), 

(24) Q(^v) = lM 2 m,V) 2 , 

where 

/(£, rj) = 1 + (1 - Co)(€ 2 m - Him) + (&<?? - 6<z 2 °) - (k + i)(i - Co) - — , 

Due to the shifting and squaring we now need to look at {Q = \^ 2 } = {K = 0}. 
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6.2. Level sets "close" to the moon or the earth are contact. Wc shall check that the 
Liouville vector field X = rid v is transverse to one connected component of the level set {Q = 
for points (£,77) with |ry|(l — £0) < e - Indeed, 

x(Q) = \v\ 2 f(t,v) 2 + \v\ 2 m,v)vdr,m,v) = 

(1-M)(1 -Co) \ 



2Q + N 2 /(e,r?) (1- 6X6% -6%) -7?^- 

V 1^7(1 - Co) + ^0 + g° - g- 1 1 y 

To estimate the last term, let us first make a few observations. First of all, let us argue that we get 
an upper bound for \t]\ from the equality Q = |/j 2 . Indeed, the term 77) I can be estimated as 



1 + (1 - £0X6771 - + (65? - 6<Z 2 °) -(k+h(l- 6) - - { ]/ ){ -} 

1 - (1 - mm m°\ + (1 - &>xi*i -\- — -v 1 ^ n — rr) 

2 |t7(1-6)+67o + <Z°-<Z 1 I 
Using |g | = 1 — fj,, |£| = 1, and (1 — 6)1^1 < e, we obtain 

1/(6 77)| > 1 - e - (1 - /x) + (1 - £ ) ( |fc| - J - — 77 J ^ f 

In the last step we have argued that 

- 1 -4 T <2(1- M ) 



> 



|77(i-6) + 67o + <7°-<z 1 l k-q 1 

is bounded on a small neighborhood of q° (we have \y\ = \q — q°\ < e). Indeed, choosing e < 1/2 is 
enough for this. In that case | <r/ — <r/ 1 1 > 1/2, so we see that the claim holds. Choose the energy level 
k appropriately. 

As Q = \ /j 2 , we obtain the promised upper bound for 

^ = <3>i|,l*£, 

so we get 

H < 2. 

Let us then continue with our earlier computation to check that the level sets are contact. We may 
write 

X(Q)> 2(3 -1771(17711/(6 77)1)77^/(^7?). 
Observe that |?7| < 2, |77||/(6t?)| = V%Q = A 1 - The remaining term can also be estimated. 



\vd v m,v)\ 



(1-6)^-6^-77^ — n 

|77(i-6) + C77o + g -9 1 l 



< e + (1 — fi)ed v — 



|7?(i-6) + ^o + g°-g 1 | 

The latter term can obviously be bounded by some constant C on a compact set away from the 
singular point. As \y\ < e, this holds true, so finally, 

X(Q) > M 2 -2/ J e(l + (l-/ J )C). 

Now we see that for e small enough we have X(Q) > on the connected component of the level set 
we are interested in. 

Remark 6.2. It is important to observe that the Liouville vector field X is actually the natural 
Liouville vector field on the cotangent bundle, i.e. rjd v . 
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Combining this with Proposition |5.1| and the computations related to that proposition, we find 
the following. 

Proposition 6.3. Let E < H(L\) and consider the component E of the level set {H = E} cor- 
responding to either q° or q . Then E can be regularized to form the closed 3-manifold E = MP 3 . 
Furthermore, E carries a natural contact structure which is isomorphic to the standard contact 
structure on RP 3 . Finally, the sublevel set {H < E} provides a filling of E, which extends to a 
strong symplectic filling o/E. The dynamics are preserved in this regularization, i.e. Reeb orbits in 
E correspond to Reeb orbits in E. 



Proof. In order to combine Proposition 5.1 with the above computations, note that the Liouville 
vector field on the regularization and the one used in Proposition |5 . 1 1 are both the natural Liouville 
vector fields on the cotangent bundle, i.e. both defined by the equation ix^can — ^can- This means 
in particular that they must coincide after the coordinate change. Hence we obtain a global Liouville 
vector field for the regularized planar circular restricted three body problem. This shows that the 
regularized level set E is contact, and diffcomorphic to MP 3 , since the complete regularized level set 
is filled by the set {Q < \p?} which can be identified with T< 1 S 2 . 

Therefore we only need to check that we obtain the standard contact structure on E = KP 3 = 
ST*S 2 . This can be done directly, but in dimension 3 we can alternatively argue as follo ws. F irst 

" On 



2.2 



of all, E is fillable, so it is tight by the Eliashberg-Gromov tightness theorem, see Theorem 
the other hand, we also know that MP 3 has a unique tight contact structure by Theorem 2.4| Hence 



the contact structure on E induced by the Liouville field X is isomorphic to the standard contact 
structure on MP 3 . □ 

7. Connected sum 

7.1. Interpolating Liouville vector fields. Let us start by describing the setup. Write q° = 
( — (1 — /i,0)) and q 1 — (/i,0) for the location of the moon and earth, respectively. Consider the 
Hamiltonian 

H{q,p) := l -\p+ Jq\ 2 + V cS (q), 



where 



and 



J 



1 

-1 



T/ 1 2 A* 1 - M 



^ \q-q°\ |<?-3T 
Following computations of Conley, [5J, we can expand the effective potential around the Lagrange 
points q L , which are the critical points of V e g. We shall write 

V cS (q) = Q(q-q L )+R(q-q L ), 

where Q is a quadratic form in q — q L given by 

1 / -2p 



Q =2 V P 

and R(q — q L ) is a rest term of higher order than (q — q L ) 2 . 

We now expand the Hamiltonian H(q,p) around the critical point (q,p) = (q L , —Jq L ). We obtain 

H(q, P ) = \\v+ Jq\ 2 + Q(q - q L ) + R(q - q L ) - \\{p + Jq L ) + J(q - q L )\ 2 + Q(q - q L ) + R(q - q L ), 

which we shall rewrite as 

H(q,p) = Q(q -q L ,p+ Jq L ) + R(q - q L ), 
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using the quadratic form Q given by 

Q = 



V 



-2p 



1 





p 

-1 





1\ 






1 / 



/ 














\ 


( d qi 


\ 







b 










dq 2 












1 -a 







dpi 




V 











1 - b 


) 


V d P2 


J 



Hence we consider the Hamiltonian 

Q(q,p)+R(q), 
for q,p near 0. We have the Weinstein-like vector field 

Ya.b = (qi,q 2 ,Pi,P2) 
Then we have the following lemma. 

Lemma 7.1. The vector field Y a ^ is Liouville. Furthermore, for a sufficiently small neighborhood 
of (q L ,-Jq L ) there exists e > 0, a < and b > such that for E € [H(L 1 )-e,H(L 1 ) + e]-{H(L 1 )} 
the vector field Y :— Y a _b is transverse to the level set H = E. 

Proof. The Liouville property can be directly checked. Furthermore, note that Y(Q) is a quadratic 
form in (q,p). For appropriate values of a and b (for instance a = — 1 and b = 1/2), we can directly 
check that Y(Q) is positive definite. In order to see this, we may use the following expression, which 
is taken from Conley [S] right below equation (4) therein, 

(25) p - 



q 



113- 



Using equation (18), we can easily check that p > 4. After that, one can verify that the eigenvalues 
of the quadratic form Y(Q) are positive. 

On the other hand, Y(R) is still of order q 2 , so for a sufficiently small neighborhood of (q L , —Jq L ) 
we can estimate this quantity by a multiple of Y(Q), 

Y(Q + R)= Y(Q) + Y(R) > Y(Q) - 1/2\Y(Q)\ > 0. 



□ 



The contact form induced by Y is given by 

ai = -a(qi - qf)dpx - bq 2 dp 2 + (1 - a)p 1 dqi + (1 - b)(p 2 - qi)dq 2 
On the other hand, by Remark 



5.6 



the vector field X = (q — q 1 )d q is Liouville and transverse to 
{(q,p)\H(q,p) — E}, for (q,p) away from the critical point L\. 



From X we obtain the induced 1-form 

«o = (q{ ~ qi)dpi - q 2 dp 2 . 

Note that 

ai - a = (1 - a) ((gi - q\)dp\ + pi<%) + (qf - ql)dpi + (1 - b) (q 2 dp 2 + (p 2 - q\ )dq 2 ) . 

is exact with primitive 

G := (1 - a){ qi - qf) Pl + (qf - q\) Pl + (1 - b)(p 2 - q{)q 2 . 

From now on we use the more convenient to use coordinates around the critical point (q L , —Jq L 
i.e. (</ ,p') = (q — q L ,p + Jq L ). With some abuse of notation (leaving out '), we obtain 

G = (1 - a)qxp\ + (qf - q\)pi + (1 - b)p 2 q 2 . 



20 



P. ALBERS, U. FRAUENFELDER, O. VAN KOERT, AND G.P. PATERNAIN 



If we denote by Xi, i — 0, 1, the Liouville vector field corresponding to the 1-form a,; and define the 
Hamiltonian vector field X G of G by 

ix G u = dG, 

then we have 

X\ = Xq + X G . 

Next we choose a cut-off function / depending on qi — j } p2 to interpolate between Xq for q\ large 
and X± for q\ close to and make the following ansatz 

X = X + XfQ- 

Remark 7.2. The dependence of / on the specific linear combination qi — ^p2 is chosen to cancel 
out some terms in a later computation. As we explain next the set {qi — ^p2 — 0} divides the 
singular energy hypersurface {H = H(Li)} into two connected components. We shall call the 
component containing q° in its closure the moon component and the component containing q 1 the 
earth component. 

Indeed, this can be seen as follows. First of all, the singular energy hypersurface {H = H(Li)} 
corresponds to {Q — 0}. We now claim that the hyperplane {qi — ^p2 = 6} intersects the quadric 
{Q = 0} in a 2-sphere that collapses to a point for 6 = 0, Indeed, inserting the equation q± — ^p2 = S 
into Q = yields 

\ ((Pi - Q2) 2 + ((P + - pS) 2 + (p- l)q 2 2 - (2p + l)ql) = 0, 

which can simplified to an equation describing a sphere, 

(Pi - Q2) 2 + (pqi ~ (p + l)Sf + (p- l)q 2 2 = (2p + l)6 2 . 

Let us now check that the Liouville field X is transverse to level sets of H close to the Lagrange 
point. By choosing a cut-off function / taking values between and 1 and vanishing outside a small 
neighborhood of the Lagrange point q L , we see that X coincides with X far away from q L . We 
checked earlier that then Xq(H) > away from the Lagrange point, see Remark 5.6 Recall also 
the definition of the Poisson bracket, 

{F,G} :=uj(X F ,X G ) = dF(X G ). 

Then we see 

X(H) = X (H)+X fG (H) = dH(X ) + {H,fG} = dH(X )-{fG,H} 
= (1 - f)dH(X ) + fdH(X ) - f{G, H} - G{f, H} 
= (1 - f)dH(X Q ) + fdH(X + X G ) + GdH(X f ). 

Note that Xq + X G = Y, so the first two terms are non-negative near the Lagrange point. For the 
last term, note that the Hamiltonian vector field of / is given by 

/did 
\ dpi p dq 2 > 



We directly compute dH(Xf). We obtain 



dH(X f )=f.[(l-% 1 + q ^ " 



p/ p \\q-q 1 \ z Iq-q 1 ^ 



Observe that the latter term vanishes in q by the equation (25) for p. As a result, we see that 
dH(Xf) ~ f'(l — -)pi, up to higher order terms. 

Now we specify the cut-off function. Choose small constants < t\ < 62- Choose a smooth 
function / such that / = 1 on [0, ei] and / = on [e2,oo), and with non-positive derivative in 
between. For later purposes we extend / symmetrically round 0. 
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Then observe that the leading order term in G is (gjp — ql)pi and the coefficient q\ — q\ is negative. 
Hence the product GdH(Xf) has a leading order term equal to 

f-(q^-ql)(l--)pl 



P> 

which is non- negative. In particular we see that the leading order terms in X(H) are all non-negative 
and at least one term is positive. Thus, they dominate the higher order terms. Hence we see that X 
is a Liouville vector field that is everywhere transverse to a component of the level set H = H(Li) + e 
provided that e is small enough. 

We have now all ingredients to prove the following lemma. 

Lemma 7.3. There exists e > and a Liouville vector field X such that X is transverse to level 
sets H = E for all E L < E < E L + e. 



Proof. As observed in Remark 7.2 the linear functional q\ — if>2, on which / depends, divides the 



level set {H = E} into a moon and earth component. 

In the above construction we already constructed a Liouville vector that is transverse to the 
earth component. Since the constructed vector field is equal to X\ on a small neighborhood of the 
Lagrange point, where / = 1, the same construction on the moon side will coincide on this ball with 
the Liouville vector field for the earth component, so we end up with a globally defined Liouville 
vector field, that is transverse to the entire level set {H = E}. □ 

7.2. Contact structures on different level sets of H. In the previous section we have shown 
that all level sets of H for energy less than H(Li) + e are of restricted contact type. On the other 



hand, Proposition 6.3 shows that level sets of H with energy less than H(Li) can be regularized 
to form a closed manifold, namely MP 3 with its natural contact structure, while preserving the 
dynamics. We pointed out in Remark |6.2| that the Liouville vector field used for the regularization 
is the natural Liouville vector field for the cotangent bundle, so it automatically coincides with the 
Liouville vector field (q — q°)d q after a suitable coordinate change. 

Now fix an energy level H(Li) < E < H(Li) + e. The above observation that the Liouville field 
is natural applies to both the earth and moon component of the level set H — E. Hence we can 
regularize the level set E = {H = E} both near q° and near q 1 in order to obtain a closed manifold 
S. 

We see directly that topologically £ = Si^X^ by observing that {H(Li) — e < H < E} provides 
a cobordism W between Si U E 2 and E. This cobordism admits H as Morse function with only one 
index 1 critical point, namely L\. Thus, we can regard E and hence E as a connected sum. 

Remark 7.4. We observe that E is symplectically tillable. The above cobordism argument shows 
that, topologically, the filling is given by the boundary connected sum T< 1 S 2 \\T< 1 S 2 . The set 
{H(Li) — e < H < E} n Bs(Li) is a symplectic handle corresponding to a 1-handle attachment as 
described in Section [2~2) A small isotopy between our setup, which is very close to the Weinstein 
model, and the actual Weinstein model induces the isotopy required to show that as contact manifolds 



E = Ei#E 2 . 



2.5 



which 



Alternatively, we can reason that, since E is tillable, it is tight. Now apply Theorem 
asserts that tight contact manifolds have a unique prime decomposition. This immediately implies 
that E = (MP 3 , £i)#(RP 3 , £2), where £1 and £2 are tight contact structures on MP 3 . Now we can 
simply use the fact that there is a unique tight contact structure on MP 3 up to isotopy just as in 
Proposition |6.3| to complete the proof of Theorem A. 



Appendix A. Computations 
Here are some detailed computations omitted in the main text of Section [5] 
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U(p, 9) = -^- 1 M - \p 2 + pcos0(l A*) 2 - 



<9C7 

= (1 — ^)psin#( r — 1 

(p 2 -2pcos6» + l)" 



89 



dU p . p-cos9 . 

— = + (1 - /i) T - p + cos6(l - p). 

dp p 2 V /y (p2_ 2 pcos0 + l)l P V W 



9 2 [/ , /sine^ 2 - + l) 2 - |(p 2 - 2/9COS0 + l)2 2psine»(p-cos0) N 

. = (1 - /i) o sin 9 

9pd6 \ ( p 2 -2pcos9 + iy } 

= (l-p) S in9( p2 - 2pCOSd + 1 - 3p{p - C ° Sd) -l) 



(p 2 -2pcos9 + l) : 
1 + p cos 9 + 1 
,(p 2 - 2pcos(9 + l)^ 



«9 2 [/ 2p (p 2 -2p cos 9 + 1)1 -^(p 2 -2p cos 9+ l)i(2p- 2 cos 9)(p- cos 9) 

dpdp ~ ~~^ + { ^~^> {pi -2p cos + 1) 3 1 

2/z p 2 - 2pcos6> + 1 - 3(p- costf) 2 ( 

P 3 (p 2 -2p cos 6> + l) I 

= H — ^(p 2 - 2pcos9 + 1 - 3p 2 + 6p cos 6 1 - 3 cos 2 6^ - 1 

p 3 (p 2 -2pcos0 + i)s v p r p ; 

2/i 1 — /x 



p 3 (p 2 -2pCOs6l+l)i 



(2p 2 -4pcos6>-l + 3cos 2 6>) -1 



d 3 U 6p (4p -4 cos 6>)(p 2 -2pcos9 + iyi 

~dpt ~ (p 2 -2pcos9 + 1) 5 

| (p 2 - 2p cos 6» + 1) i (2p - 2 cos (9) (2p 2 - 4p cos 6» - 1 + 3 cos 2 6) 
+ (1 ~ M) (p 2 -2pcos6i + l) 5 

6p / , N 4(p 2 -2 / 9Cos6»+l)-5(2p 2 -4pcos6»-l + 3cos 2 6') 

= — - I 1 - p)(p-cos6») — — T 

P 4 (p 2 - 2pcos9 + 1)5 

6p ,„ w ^-6p 2 + 12p cos 61 + 9-15 cos 2 9 

= -j-( l -v)(p-cos9) — — 

P (p 2 - 2pcos9 + 1)2 

6p „, w n . 3 + Apcos9 - 2p 2 - 5cos 2 6> 

p (p 2 — 2pcos9 + 1)2 
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U(p,0)-U(p,w) = -t^—+P+t^—+P 



1-p l+p 

-l-p+l-p + 2p(l-p)(l + p) 



(l-p)(l + p) 

— On -I- On 

= (i-m)- 



2p + 2p- 2p 3 



(l-p)(l + p) 
2(1 -P)P 3 



(l-p)(l + p) 

< 0. 



= (1- M)(2-7r^-7TT^2 



= (i-m) 



(1-p) 2 + 

2(1-P 2 ) 2 -(1 + P) 2 -(1-P) 2 

(l-p)2(l + p)2 

2 - 4p 2 + 2p 4 - 1 - 2p - p 2 - 1 + 2p - p 2 



2(1 -p) 



(l-p)2(l + p)S 

P 4 -3p 2 



(l-p) 2 (l + p) 2 



_ 2p 2 - 4pcos6< - 1 + 3cos 2 6> 
' 9( ~ (p 2 -2pcos0 + l)§ 



k(6>) = p 2 - 2pcos0 + 1. 
(4psin6< - 6cos6 , sin6>)^ - §2 ( osin6ki (2p 2 - 4pcos6* - 1 + 3 cos 2 0) 



sinfl 
sinfl 

K2 



sin # / o 2 2 

— ^ — 6p + 6p cos — 3p cos # + 9p — 6cosi 

_ ^ cos 2 6- + 2(1 - p 2 ) cos 6- + p(2p 2 - 3) 



(4p-6cos0)(p 2 -2pcos0 + l) - 5p(2p 2 -4pcos0- l + 3cos 2 9) 
4p 3 - 8p 2 cos + 4p - 6p 2 cos B + 12 p cos 2 6> - 6 cos 6 
-10p 3 + 20p 2 cos 6» + 5/?- 15/9 cos 2 6» 



6> e M. 
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2(p 2 - 1) ± y/4(l - p 2 ) 2 _ 4 p 2 (2p 2 ~) 

COS a = 

2p 

p 2 - 1 ± y/l - 2p 2 + p 4 - 2p 4 + 3p 2 
P 

p 2 - 1 ± y/-p 4 + p 2 + 1 
P 

p 2 - 1 + y/-p 4 + p 2 + 1 
P 



2p 2 - Ap 2 + 4 - 4^-p 4 + p 2 + 1 - 1 + ^ ((p 2 - l) 2 + 2(p 2 - 1)V-P 4 + P 2 + 1- P 4 + P 2 + l) 

(p 2 - 2p 2 + 2 - 2^-p 4 + p 2 + 1 + 1) 1 
-2p 4 + 3p 2 - 4p 2 y/-p 4 + p 2 + 1 + 3p 4 - 6p 2 + 3 + 6(p 2 - l)y/-p 4 + p 2 + 1 - 3p 4 + 3p 2 + 3 

(3 - p 2 - 2 V /-P 4 + P 2 + 1) I P 2 
6 - 2p 4 - (6 - 2p 2 )y/-p 4 + p 2 + 1 
(3 - p 2 - 2 v /-p 4 + p 2 + l) S p 2 
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3 + 4/3 cos (9 -2p 2 -5cos 2 6> = 3 + 4p 2 - 4 + 4^-p 4 + p 2 + l-2p 2 

P 

2p 2 -l + A^f-p 4 + P 2 + 1 



:> ' 4 2p 2 + 1 + 2(p 2 - 1)V-P 4 + P 2 + 1 - P 4 + P 2 + 1 



P 2 



-^-p 2 + 2 + 2(p 2 -l)v/-p 4 + P 2 + lj 

2p 2 - 1 + 4^-p 4 + p 2 + 1 + 5 - ^ 

P 

-10V-P 4 + p 2 + 1 + ^ V-p 4 + P 2 + 1 
P 

2p 2 + 4 - 6V-P 4 + P 2 + 1 - ^ + ^ V-P 4 + P 2 + 1 

p z p A 

2p 4 + 4p 2 - 10 + (10 - 6p 2 )v/-p 4 + P 2 + 1 
P 2 

(10 - 6p 2 ) 2 (-p 4 + p 2 + 1) - (2p 4 + 4p 2 - 10) 2 
p 2 (l0 - 2p 4 - 4p 2 + (10 - 6p 2 )\/-p 4 + p 2 + l) 

(36p 4 - 120p 2 + 100) (-p 4 + p 2 + 1) 
p 2 (l0 - 2p 4 - 4p 2 + (10 - 6p 2 )v/-p 4 + p 2 + l) 
(4p 8 + 8p 6 - 20p 4 + 8p 6 + 16p 4 - 40p 2 - 20p 4 - 40p 2 + 100) 
p 2 (lO - 2p 4 - 4p 2 + (10 - 6p 2 )\/-p 4 + p 2 + l) 
-36p 8 + 120p 6 - 100p 4 + 36p 6 - 120p 4 + 100p 2 + 36p 4 - 120p 2 + 100 
p 2 (lO - 2p 4 - 4p 2 + (10 - 6p 2 )v/-p 4 + p 2 + l) 
-4p 8 - 16p 6 + 24p 4 + 80p 2 - 100 
p 2 (lO - 2p 4 - 4p 2 + (10 - 6p 2 )y/-p 4 + p 2 + l) 
-40p 8 + 140p 6 - 160p 4 + 60p 2 
p 2 (lO - 2p 4 - 4p 2 + (10 - 6p 2 )v/-p 4 + p 2 + l) 

20(3-8p 2 + 7p 4 -2p 6 ) 
10 - 2p 4 - 4p 2 + (10 - 6p 2 )v/-p 4 + p 2 + 1 



dU _ rp (l-r)(l-p) 
dr 



3 (r-l + p). 
(r 2 + s 2 Y ({l-r) 2 + s 2 Y 



dU sp s(l — p) 

Os ( r 2 + s 2)l ((1 _ r )2 + s 2) I 
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dU dU (r 2 + .S 2 )p (I-/*) , . 2 , 

r 7T +s lT = ^ + 5-(r(r-l) + s 2 )+r(l-r-Ai) 

M r + 3 (r 2 + s 2 -r)+r(l-r-n)-s 2 



s 2 



(r 2 + s 2 ) 5 ((l-r) 2 + s 2 ) 

(Pi +s) 2 + (p 2 + l-r- p) 2 = 2{k - U). 

A* 1 — A* / 2 2 \ 

— === -p 2 r+pisH r (r +s - r) 

Vr 2 + s 2 ((l-r) 2 + s 2 ) 5 

= , ; - (P2 + 1 - r - p;)r + r(l - r - p) + (p 1 + s)s - s 2 
vr + s l 

+ ^ Hr 2 + s 2 -r) 

((l-r) 2 + s 2 ) 5 

> -4==+r(l-r-v)- S 2 + - N 3 (r 2 + . 2 -r) 

Vr 2 + s 2 ((l-r) 2 + s 2 ) 5 

-V 7 ? -2 + s 2 ^/(pi + s) 2 + (p 2 + 1 - r - a*) 2 

= r ¥+ s ?- Vr 2 + aV2(fe-t/)- 
or as 

p 5 - (3 - p)p 4 + (3 - 2p)p 3 - pp 2 +2np-fi = 0. 

p 5 - 3p 4 + 3p 3 
M ~ p 4 - 2p 3 - p 2 + 2p - 1 ' 

2p 3 (p 2 -3p + 3) 
(p 4 -2p 3 -p 2 + 2p-l)p 3 

(p 4 - 2p 3 - p 2 + 2p - 1 + p 5 - 3p 4 + 3p 3 ) (6 - 2p 4 - (6 - 2p 2 ) ^-p 4 + p 2 + l) 
(p 4 - 2p 3 - p 2 + 2p - 1) (3 - p 2 - 2 V /-P 4 + P 2 + 1) § P 2 
2p 2 (p 2 - 3p + 3) (3 - p 2 - 2y/-p 4 + p 2 + 1) 1 
(p 4 - 2p 3 - p 2 + 2p - 1) (3 - p 2 - 2 v /-p 4 + p 2 + l) i p 2 



(p 5 - 2p 4 + p 3 - p 2 + 2p - 1) (6 - 2p 4 - (6 - 2p 2 ) v /-p 4 + (0 2 + i) 



(p 4 - 2p 3 - p 2 + 2p - 1) (3 - p 2 - 2^-p 4 + p 2 + 1) 1 p 2 

A.l. Proof that W(po,0o) is negative. Observe that 

d 5 - 2d 4 + d 3 - d 2 + 2d - 1 = (d 2 + d + l)(d - l) 3 , 

and that 

d 4 - 2d 3 - d 2 + 2d - 1 = (d - 1/2) 4 - 5/2(d - 1/2) 2 - 7/16. 

The latter can be checked to be negative on the interval [0, 1] We now rewrite everything to see that 
this expression is negative. The main trick which we shall repeat over and over again, is the simple 
identity 

a - Vb = — p. 

a + Vb 
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For instance, we see that 

j2\2 



3 - d 2 + 2\/l + d 2 - d 4 3 - d 2 + 2\/l + d 2 - d 4 

Hence we obtain the following: 

2d 2 (d 2 - 3d + 3) (3 - d 2 - 2V-d 4 + d 2 + l) 1 



(d 4 - 2d 3 - d 2 + 2d - 1) (3 - d 2 - 2V-d 4 + d 2 + l)*d 2 
(d 5 - 2d 4 + d 3 - d 2 + 2d - 1) (6 - 2d 4 - (6 - 2d 2 )V-d 4 + d 2 + l) 
(d 4 - 2d 3 - d 2 + 2d - 1) (3 - d 2 - 2yJ-d 4 + d 2 + l)h 2 

2(d 2 1 31 ( 5(i-d) 2 (i+d) a y 



5/2 



5/2 



((d - 1/2) 4 - 5/2(d - 1/2) 2 - 7/16) ( 3 XZi + T d < y 
(d 2 + d + l)(d - 1)3 (-6 1 ;; 1 ^ 4 + 2V1 + d 2 - d 4 - 2d 2 ) 
((d - 1/2) 4 - 5/2(d - 1/2) 2 - 7/16) (^BSf) 
Note here that d e [0, 1], so we can pull out (1 — d) 5 in the expression above. We obtain 

2(d 2 - 3d + 3)(1 - df(l + df (3 _ j+ Xp=?)v» 
((d - 1/2) 4 - 5/2(d - 1/2) 2 - 7/16) (1 - d) B (l + ^ V^T*-*)'/* 

2(d 2 + d + 1)(1 - df 1 ;; 1 ^ 4 (-2 + VI + d 2 - d 4 + d 2 ) 
+ ((d - 1/2) 4 - 5/2(d - 1/2) 2 - 7/16) (1 - d)«(l + «9 V* + Xp-*)»/» 



= 2 



(d 2 -3d + 3)(l-d) 5 (l + d) E 



(3-d 2 +2Vl+d 2 -<i 4 ) 5/2 



((d - 1/2) 4 - 5/2(d - 1/2) 2 - 7/16) (1 - d) B (l + ^ V^T*-*)'/* 
(^ +d+ i)(i- rf )^(i + ^(-i +i+Vi ^ 2 _ d4 



, l + Vl+d 2 -d 4 

((d - 1/2) 4 - 5/2(d - 1/2) 2 - 7/16) (1 - d)«(l + d) V* + Xp-*)»/» 

C5/2/1 i j\3 d 2 -3d+3 , (rf 2 +d+l)d 2 (rf 2 +rf+l) 

2 _ l 1 ^"! (3-d 2 +27l + d 23 rf 3 ') !s/2 (l+Vl+rf 2 -d 4 ) 2 l + Vl+d 2 -d 4 

((d - 1/2) 4 - 5/2(d - 1/2) 2 - 7/16) (1 + d) 3 (3 _ JJ+ XV*)^ 

Now observe that the denominator is negative for d G [0, 1]; there are no zeroes for d € [0, 1] and 
the function (d — 1/2) 4 — 5/2(d — 1/2) 2 — 7/16 is clearly negative. Furthermore, the denominator is 
positive for the following reason. We shall estimate the sum of the first and last term, 

5 5 / 2 (l + d) 3 (d 2 -3d + 3) (d 2 +d+l) 



(3 - d 2 + 2y/l + d 2 - d 4 ) 5 /2 i + y/l + d 2 - d 4 

5^ 1 

. (3 - d 2 + 2y/l + d 2 - d 4 ) 5 / 2 1 + VI + d 2 - d 4 



(d 2 + d + 1) 



We see that this estimate holds, since for d € [0, 1] we have d 2 +d+ 1 < (1 + d) 3 , and d 2 — 3d + 3 > 1; 
indeed, (d 2 — 3d + 3)' = 2d — 3 is negative on [0, 1], so d 2 — 3d + 3 is minimal in d = 1. Finally 
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observe that j-^==== < 1, whereas 

55/2 



(3 - d 2 + 2V1 + d 2 - d 4 f/ 2 
The latter can be seen by observing that (3 — d 2 + 2\f\ + d 2 — d 4 ) has derivative 

-2d (l - /~ 2rf2 ) < -2d(- ) < 0, 

so 

5 5/2 



(3 - d 2 + 2y/l + d 2 - d 4 ) 5 / 2 
is minimal for d = 0, where its value is 1. This shows the claim. 
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